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Abstract 

In this paper we extend the existing literature on the asymptotic behaviour 
of the partial sums and the sample covariances of long memory stochastic 
volatility models in the case of infinite variance. We also consider models 
with leverage, for which our results are entirely new in the infinite variance 
case. Depending on the interplay between the tail behaviour and the intensity 
of dependence, two types of convergence rates and limiting distributions can 
arise. In particular, we show that the asymptotic behaviour of partial sums is 
the same for both LMSV and models with leverage, whereas there is a crucial 
difference when sample covariances are considered. 
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1. Introduction 

One of the standardized features of financial data is that returns are uncorrelated, but their 
squares, or absolute values, are (highly) correlated, a property referred to as long memory (which 
will be later defined precisely) . A second commonly accepted feature is that log-returns are heavy 
tailed, in the sense that some moment of the log-returns is infinite. The last one we want to 
mention is leverage. In the financial time series context, leverage is understood to mean negative 
dependence between previous returns and future volatility (i.e. a large negative return will be 
followed by a high volatility). Motivated by these empirical findings, one of the common modeling 
approaches is to represent log-returns {Y{\ as a stochastic volatility sequence Yj, = ZiOi where {Zi} 
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is an i.i.d. sequence and {erf} is the conditional variance or more generally a certain process which 
stands as a proxy for the volatility. In such a process, long memory can only be modeled through 
the sequence {o~i}, and the tails can be modeled either through the sequence {Zi} or through 
{<7i}, or both. The well known GARCH processes belong to this class of models. The volatility 
sequence {o~i} is heavy tailed, unless the distribution of Zq has finite support, and leverage can 
be present. But long memory in squares cannot be modeled by GARCH process. The FIGARCH 
process was introduced by [3] to this purpose, but it is not known if it really has a long memory 
property, see e.g. [15]. 

To model long memory in squares, the so-called Long Memory in Stochastic Volatility (LMSV) 
process was introduced in [7], generalizing earlier short memory version of this model. In this 
model, the sequences {Zi} and {o~i} are fully independent, and {(Ji} is the exponential of a Gaussian 
long memory process. Tails and long memory are easily modeled in this way, but leverage is absent. 
Throughout the paper, we will refer to this process as LMSV, even though we do not rule out the 
short memory case. 

In order to model leverage, [26] introduced the EGARCH model (where E stands for exponen- 
tial), later extended by [6] to the FIEGARCH model (where FI stands for fractionally integrated) 
in order to model also long memory. In these models, {Zi} is a Gaussian white noise, and {o~i\ is 
the exponential of a linear process with respect to a function of the Gaussian sequence {Zi}. [32] 
extended the type of dependence between the sequences {Zi} and {Xi} and relaxed the Gaussian 
assumption for both sequences, but assumed finite moments of all order. Thus long memory and 
leverage are possibly present in these models, but heavy tails are excluded. 

A quantity of other models have been introduced, e.g. models of Robinson and Zaffaroni [29], 
[30] and their further extensions in [28]; LARCH(oo) processes [19] and their bilinear extensions 
[20], and LARCH + (oo) [31]; to mention a few. All of these models have long memory and some 
have leverage and allow for heavy tails. The theory for these models is usually extremely involved, 
and only the asymptotic properties of partial sums are known in certain cases. We will not consider 
these models here. In [18] the leverage effect and long memory property of a LARCH(oo) model 
was studied thoroughly. 

The theoretical effect of long memory is that the covariance of absolute powers of the returns 
{Yi} is slowly decaying and non summable. This induces non standard limit theorems, such as 
convergence of the partial sum process to the fractional Brownian motion or finite variance non 
Gaussian processes or even Levy processes. In practice, long memory is often evidenced by sample 
covariance plots, showing an apparent slow decay of the covariance function. Therefore, it is of 
interest to investigate the asymptotic behaviour of the sample mean or of the partial sum process, 
and of the sample variance and covariances. 

In the case where <7i = o~{Xi), {Xi} is a stationary Gaussian process with summable covariances 
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and <j{x) = exp(cc), the asymptotic theory for sample mean of LMSV processes with infinite 
variance is a straightforward consequence of a point process convergence result in [14]. The 
limit is a Levy stable process. [32] considered the convergence of the partial sum process of 
absolute powers of generalized EGARCH processes with finite moments of all orders and showed 
convergence to the fractional Brownian motion. To the best of our knowledge, the partial sum 
process of absolute powers has never been studied in the context of heavy tails and long memory 
and possible leverage, for a general function a. 

The asymptotic theory for sample covariances of weakly dependent stationary processes with 
finite moments dates back to Anderson, see [1]. The case of linear processes with regularly varying 
innovations was studied in [10] and [11], for infinite variance innovation and for innovations 
with finite variance but infinite fourth moment, respectively. The limiting distribution of the 
sample covariances (suitably centered and normalized) is then a stable law. These results were 
obtained under conditions that rule out long memory. For infinite variance innovation with tail 
index a 6 (1,2), these results were extended to long memory linear processes by [24]. The 
limiting distributions of the sample covariances are again stable laws. However, if a € (2,4), [21] 
showed that as for partial sums, a dichotomy appears: the limiting distribution and the rate of 
convergence depend on an interplay between a memory parameter and the tail index a. The limit 
is either stable (as in the weakly dependent or i.i.d. case) or, if the memory is strong enough, the 
limiting distribution is non Gaussian but with finite variance (the so-called Hermitc-Roscnblatt 
distributions). If the fourth moment is finite, then the dichotomy is between Gaussian or finite 
variance non Gaussian distributions (again of Hermite-Rosenblatt type); see [22], [21, Theorem 
3.3] and [34]. 

The asymptotic properties of sample autocovariances of GARCH processes have been studied 
by [4]. Stable limits arise as soon as the marginal distribution has an infinite fourth moment. 
[14] studied the sample covariance of a zero mean stochastic volatility process, under implicit 
conditions that rule out long memory, and also found stable limits. [25] (generalized by [23]) 
studied partial sums and sample variance of a possibly nonzero mean stochastic volatility process 
with infinite variance and where the volatility is a Gaussian long memory process (in which case 
it is not positive but this is not important for the theoretical results). They obtained a dichotomy 
between stable and finite variance non Gaussian limits, and also the surprising result that when 
the sample mean has a long memory type limit, then the studentized sample mean converges in 
probability to zero. 

The first aim of this article is to study asymptotic properties of partial sums, sample variance 
and covariances of stochastic volatility processes where the volatility is an arbitrary function of a 
Gaussian, possibly long memory process {X{\ independent of the sequence {^}, which is a heavy 
tailed i.i.d. sequence. We refer to these processes as LMSV processes. The interest of considering 
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other functions than the exponential function is that it allows to have other distributions than the 
log- normal for the volatility, while keeping the convenience of Gaussian processes, without which 
dealing with long memory processes becomes rapidly extremely involved or even intractable. The 
results we obtain extend in various aspects all the previous literature in this domain. 

Another important aim of the paper is to consider models with possible leverage. To do this, we 
need to give precise assumptions on the nature of the dependence between the sequences {Zi} and 
{X;}, and since they are related in the process {Y{\ through the function a, these assumptions 
also involve the function a. We have not looked for the widest generality, but the functions a 
that we consider include the exponential functions and all symmetric polynomials with positive 
coefficients. This is not a severe restriction since the function a must be nonnegative. Whereas 
the asymptotic theory for the partial sums is entirely similar to the case of LMSV process without 
leverage, asymptotic properties of sample autocovariances may be very different in the presence 
of leverage. Due to the dependence between the two sequences, the rates of convergence and 
asymptotic distribution may be entirely different when not stable. 

The article is organized as follows. In Section 2 we formulate proper assumptions, as well as 
prove some preliminary results on the marginal and multivariate tail behaviour of the sequence 
{Yi}. In Section 3, we establish the limit theory for a point process based on the rescaled sequence 
{Yi}. This methodology was first used in this context by [14] and our proofs are closely related 
to those in this reference. Section 4 applies these results to obtain the functional asymptotic 
behaviour of the partial sum process of the sequences {Yi} and of powers. In Section 5 the 
limiting behaviour of the sample covariances and autocorrelation of the process {Yi} and of its 
powers is investigated. Proofs are given in Section 6. In the Appendix we recall some results on 
multivariate Gaussian processes with long memory. 

A note on the terminology 

We consider in this paper sequences {Yi} which can be expressed as Yi = Zia(Xi) = Zi<Ji, 
where {Zi} is an i.i.d. sequence and Zi is independent of Xi for each i. Originally, SV and LMSV 
processes refer to processes where the sequences {Zi} and {<t,} are fully independent, at = a(Xi), 
{Xi} is a Gaussian process and a(x) = exp(x); see e.g. [7], [8], [14]. The names EGARCH 
and FIEGARCH, introduced respectively by [26] and [6], refer to the case where cr(x) — cxp(a;) 
and where {Xi} is a non Gaussian process which admits a linear representation with respect 
to an instantaneous function of the Gaussian i.i.d. sequence {Zi}, with dependence between the 
sequences {Zi} and {X;}. [32] still consider the case <j(x) = exp(a;), but relax the assumptions on 
{Zi} and {Xi}, and retain the name EGARCH. The LMSV processes can be seen as border cases 
of EGARCH type processes, where the dependence between the sequences {Zi} and {Xi} vanishes. 
In this article, we consider both LMSV models, and models with leverage which generalize the 
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EGARCH models as defined by [32]. In order to refer to the latter models, we have chosen 
not to use the acronym EGARCH or FIEGARCH, since these models were defined with very 
precise specifications and this could create some confusion, nor to create a new one such as 
GEGARCH (with G standing twice for generalized, which seems a bit too much) or (IV)LMSVwL 
(for (possibly) Infinite Variance Long Memory Stochastic Volatility with Leverage). Considering 
that the main feature which distinguishes these two classes of models is the presence or absence 
of leverage, we decided to refer to LMSV models when leverage is excluded, and to models with 
leverage when we include the possibility thereof. 



2. Model description, assumptions and tail behaviour 

Let {Zi,i £ Z} be an i.i.d. sequence whose marginal distribution has regularly varying tails: 

Hm F{Z °^ =P, Hm r ^<-'K l-p, (1) 

x^+oo X - a L(x) x-y+oo x~ a L(x) 

where a > 0, L is slowly varying at infinity, and /3 £ [0, 1]. Condition (1) is referred to as the 
Balanced Tail Condition. It is equivalent to assuming that P(|Zo| > x) = x~ a L{x) and 

x^ + oo F(\Zq\ > x) x^+oo ¥(\Z \ > x) 

We will say that two random variables Y and Z are right-tail equivalent if there exists c £ (0, oo) 
such that 

V(Y > x) 
hm — - = c . 

x^+oo F(Z > x) 

If one of the random variables has a regularly varying right tail, then so has the other, with the 
same tail index. The converse is false, i.e. two random variables can have the same tail index 
without being tail equivalent. Two random variables Y and Z are said to be left-tail equivalent if 
—Y and — Z arc right-tail equivalent, and they are said to be tail equivalent if they are both lcft- 
and right-tail equivalent. 

Under (1), if moreover E[|Zo| Q ] = oo, then Z\Z^ is regularly varying and (see e.g. [11, 
Equation (3.5)]) 

lim y Z ° >X \=Q, 
x^+oo ¥{Z a Zi > x) 

lim ^f^W + (l-/3) 2 . 

a^+oo F(\ZiZ 2 \ > X) H V H> 

For example, if (1) holds and the tail of \Z \ has Pareto-type tails, i.e. P(|Z | > x) ~ cx~ a as 
x — > +oo for some c > 0, then E [|Zo| Q ] = oo. We will further assume that {Xi} is a stationary 
zero mean unit variance Gaussian process which admits a linear representation with respect to an 
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i.i.d. Gaussian white noise {rji} with zero mean and unit variance, i.e. 

oo 

Xi = c Mi-i ( 2 ) 

3=1 

with J^jLi cj = 1- We assume that the process {Xi} either has short memory, in the sense 
that its covariancc function is absolutely summable, or exhibits long memory with Hurst index 
H € (1/2, 1), i.e. its covariance function {p n } satisfies 

oo 

p n = cov(X ,X n ) = CjC j+n = n 2H ' 2 i{n) , (3) 
3=1 

where I is a slowly varying function. 

Let a be a deterministic, nonnegative and continuous function defined on R. Define cr,; = cr(Xj) 
and the stochastic volatility process {Yi} by 

Yi = = a{Xi)Zi . (4) 

At this moment we do not assume independence of {77^} and {Zi}. Two special cases which we 
are going to deal with are: 

• Long Memory Stochastic Volatility (LMSV) model: where {77^} and {Z{\ arc independent. 

• Model with leverage: where {(j]i,Zi)} is a sequence of i.i.d. random vectors. For fixed i, Zi 
and Xi are independent, but Xi may not be independent of the past {Zj, j < i}. 

Both cases are encompassed in the following assumption which will be in force throughout the 
paper. 

Assumption 1. The Stochastic Volatility process {Yi} is defined by 

Yi = GiZi , 

where er,; = a(Xi), {Xi} is a Gaussian linear process with respect to the i.i.d. sequence {r/i} of 
standard Gaussian random variables such that (2) holds, a is a nonnegative function such that 
P(cr(a?7o) > 0) = 1 for all a 7^ 0, {{Zi,rji)} is an i.i.d. sequence and Zq satisfies the Balanced Tail 
Condition (1) with E[|Z | Q ] = 00. 

Let Ti be the sigma-field generated by rjj, Zj, j < i. Then the following properties hold. 

• Zi is ^-measurable and independent of J^i—i; 

• Xi and (Tj are J-"i_i-measurablc. 

We will also impose the following condition on the continuous function a. There exists q > 
such that 

sup E[<j q (jXo)} < 00 . (5) 

0<7<1 
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It is clearly fulfilled for all q, q' if a is a polynomial or u{x) = exp(x) and Xq is a standard Gaussian 
random variable. Note that if (5) holds for some q > 0, then, for q' < q/2, it holds that 



sup E 

0<7<1 



< oo , s = 1, 2, 



2.1. Marginal tail behaviour 

If (5) holds, then clearly E[ct 9 (Xo)] < oo. If moreover q > a, since Xi and Zi are independent 
for fixed i, Breiman's Lemma (see e.g. [27, Proposition 7.5]) yields that the distribution of Y$ is 
regularly varying and 

I- T° > X \ = hm ^ < ~ X \ = E K (X o)] . (6) 
x^+oc F(Z > x) x^ +oa ¥(Z <-x) 1 v 7J v 

Thus we see that there is no effect of leverage on marginal tails. Define 

a n = M{x:P{\Y \>x) <l/n} . (7) 

Then the sequence a n is regularly varying at infinity with index l/a. Moreover, since a is 
nonncgative, Zq and Yq have the same skewness, i.e. 

lim nF(Y > a n ) = 1 - lim n¥(Y < -a n ) = P ■ 

2.2. Joint exceedances 

One of the properties of heavy tailed stochastic volatility models is that large values do not 
cluster. Mathematically, for all h > 0, 

P(|F | > x, \Y h \ >x) = o(P(\Y \ > x)) . (8) 

For the LMSV model, conditioning on oo, ah yields 

P(ra>xini>x ) = 

:r-> + oo P 2 (|Z | > XJ 

if (5) holds for some q > 2a. Property (8) still holds when leverage is present. Indeed, let Fz 
denote the distribution function of Zq and Fz = 1 — Fz- Recall that Th-\ is the sigma-field 
generated by rjj, Zj, j < h — 1. Thus, Yq and Xh are measurable with respect to J~h-i, and Zh is 
independent of IFh-i- Conditioning on Fh-i yields 

P(Yo >x,Y h >x) = E[F z (x/a h )l {Yo>x} ] . 

Next, fix some e > 0. Applying Lemma 6.2, there exists a constant C such that for all x > 1, 

~F Z ( X /(T h ). 



¥(Y Q >x,Y h >x) _ E 
P(Z > x) 



-l{r >z} 



< CE [(1 V ( r h ) Q+e l {yo>:c} ] 



If (5) holds for some q > a, and e is chosen small enough so that a + e < then by bounded 
convergence, the latter expression is finite and converges to as x — > +oo. 
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2.3. Products 

For the LMSV model, another application of Breiman's Lemma yields that Y^Yh is regularly 
varying for all h. If (5) holds for some q > 2a, then 

iH? P^y 7 ^ A = E foOfc) , lim p7y y ~ r = MP*°h) . 10 

x^+oo F(ZqZi > XJ z->+oo T\ZqZ\ < —X) 

For further reference, we gather in a Lemma some properties of the products in the LMSV case, 
some of which are mentioned in [14] in the case cr(x) = exp(x). 

Lemma 2.1. Let Assumption 1 hold and let the sequences {rji} and {Z{\ be mutually independent. 
Assume that (5) holds with q > 2a. Then YqY\ is tail equivalent to ZqZ\ and has regularly varying 
and balanced tails with index a. Moreover, for all h > 1, there exist real numbers d+(h), d-{h) 
such that 

lim nYoY h >x) h y<-»U W , 
x-+ooP(\Y Y L \>x) +V ' *-><» P(|F Fi| > x) y ' y 1 

Let b n be defined by 

b n = inf{a: : P(|y n| > x) < l/n] . (12) 
The sequence {b n } is regularly varying with index l/a and 

a n = o(b n ) ■ (13) 

For all i =/= j > 0, it holds that 

lim n¥(\Y \ > a n x , \Y Q YA > b n x) = , (14) 

n— too 

lim nF(\Y Yi\ > b n x , \Y Y,\ > b n x) = . (15) 

n—too 

The quantities d+{h) and d-(h) can be easily computed in the LMSV case. 

o2 . /-i a , 2 ,E[a a (X )a a (X h )} , n ,_ oan Q ,E[a a (X )a a (X h )] 



d + {h) = {(3' + (1 - „ " „ " , <*_(/») = 2/3(1 - /?) 



When leverage is present, many different situations can occur, obviously depending on the 
type of dependence between Zq and 770, and also on the function a. We consider the exponential 
function o~(x) = exp(a;), and a class of subadditive functions. In each case we give an assumption 
on the type of dependence between Zq and 770 that will allow to prove our results. Examples are 
given after the Lemmas. 

Lemma 2.2. Assume that o~{x) = cxp(a;) and cxp(krjo) Zq is tail equivalent to Zq for all k € K. 
Then all the conclusions of Lemma 2.1 hold. 
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Lemma 2.3. Assume that the function a is subadditive, i.e. there exists a constant C > such 
that for all x,y £ R, a(x + y) < C{a(x) + o~(y)}. Assume that for any a, b > 0, cr(a£ + br/o)Zo 
is tail equivalent to Zq, where £ is a standard Gaussian random variable independent 0/770, and 
cr(br)o)Zo is either tail equivalent to Zq or M[{<7(brjo)\Zo\} q ] < 00 for some q > a. Then all the 
conclusions of Lemma 2.1 hold. 

Example 1. Assume that Zq = I770 1 1 / a! fTo with a > 0, where Uq is independent of 770 and 
E[|[/o| 9 ] < 00 for some q > a. Then Zq is regularly varying with index —a. 

• Case <t(x) = exp(x). For each c > 0, Zoexp(c7yo) is tail equivalent to Zq. See Lemma 6.1 
for a proof of this fact. 

• Case cr(x) = x 2 . Let q' e (a, q A {a/ (1 - 2a) + }). Then 

E[cr"'(b Vo )\Z \"'} = b^'n^'^I^U^'] < 00 . 
Furthermore, let £ be a standard Gaussian random variable independent of 770 and Zq. Then, 
a« + M^o = a 2 eZQ + 2a^sign(77o)|7, o r- 1/Q [/o + b 2 \ m \ 2 - l ' a UQ . 

Since £ is independent of Zq and Gaussian, by Breiman's lemma, the first term on the right- 
hand side of the previous equation is tail equivalent to Zq. The last two terms have finite 
moments of order q' for some q' > a and do not contribute to the tail. Thus the assumptions 
of Lemma 2.3 are satisfied. 

Example 2. Let Zq have regularly varying balanced tails with index —a, independent of 770. Let 
^i(-) and 4 , 2(') be polynomials and define Zq = 1(770) + ^2(^0)- Then, by Breiman's Lemma, 
Zq is tail equivalent to Zq, and it is easily checked that the assumptions of Lemma 2.2 are satisfied 
and the assumptions of Lemma 2.3 are satisfied with a being any symmetric polynomial with 
positive coefficients. We omit the details. 

3. Point process convergence 

For s = 0, . . . , h, define a Radon measure A s on [—00, co] \ {0} by 

Ao(dz) = a {[Sx-^l^ix) + (1 - j 8)(-i)-° , - 1 l ( _ 00 ,o)(a:)} dx , 
X s {dx) = a {d+(s)x^ a ^ 1 l { Q t0c) {x) + d_(s)(-a;) _Q_1 l ( „ OO:0) (a:)} dx , 

where d±(s) are defined in (11). For s = 0,...,h, define the Radon measure v s on [0,1] x 
[-oo,oo]\{0} by 

v s {At,dx) =dtA s (da;) . 
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Set Y U:i = (a^Yi, b^YiYi+i, . . . , b n l Y i Y i+ h) 1 where a n and b n are defined in (7) and (12) 
respectively, and let N n be the point process defined on [0, 1] x ([— oo, oo]' l+1 \ {0}) by 

n 

N n = ^ 6(i/n,Y„ ti ) i 
i=l 

where 5 X denotes the Dirac measure at x. 

Our first result is that for the usual univariate point process of exceedances, there is no effect 
of leverage. This is a consequence of the asymptotic independence (8). 

Proposition 3.1. Let Assumption 1 hold and assume that a is a continuous function such that 
(5) holds with q > a. Then YH=i${i/n,Yi/a n ) converges weakly to a Poisson point process with 
mean measure v$. 

For the multivariate point process N n , we consider first LMSV models and then models with 
leverage. 

3.1. Point process convergence: LMSV case 

Proposition 3.2. Let Assumption 1 hold and assume that the sequences {r/i} and {Zi} are 
independent. Assume that the continuous volatility function a satisfies (5) for some q > 2a. 
Then 

h oo 
i=0 k=l 

where YlkLi ^(t k ,jk o)' ■ • • > ^2k=i h) are independent Poisson processes with mean measures 

uq, ... , Vh, and e, € M. h+1 is the i-th basis component. Here, denotes convergence in distribution 
in the space of Radon point measures on (0, 1] x [— oo, oo]' l+1 \{0} equipped with the vague topology. 

3.2. Point process convergence: case of leverage 

Proposition 3.3. Let Assumption 1 hold. Assume that a(x) = exp(x) and Zoexp(c?yo) is tail 
equivalent to Zq for all c. Then the convergence (16) holds. 

Proposition 3.4. Let Assumption 1 hold. Assume that the distribution of (Zq,t]q) and the 
function a satisfy the assumptions of Lemma 2.3 and moreover 

\a{x + y)~a(x + z)\< C(a(x) V l){(a(y) V 1) + (a(z) V l)}\y - z\ . (17) 

Assume that condition (5) holds for some q > 2a. Then the convergence (16) holds. 

The condition (17) is an ad- hoc condition which is needed for a truncation argument used in the 
proof. It is satisfied by all symmetric polynomials with positive coefficients. (The proof would not 
be simplified by considering polynomials rather than functions satisfying this assumption.) 
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4. Partial Sums 

Define 

[nt] [nt] 

s n (t) = Y,Yi, swt) = £w- 

i=l i=l 

For any function g such that E[g 2 (?7o)] < oo and any integer q > 1, define 

J 9 ( 5 )=E[F 9 (%)5(%)] , 

where iJ g is the g-th Hermite polynomial. The Hermite rank r(g) of the function g is the smallest 
positive integer r such that J T (g) 7^ 0. Let R Tj h be the so-called Hermite process of order r 
with self-similarity index 1 — t(1 — iJ). See [2] or Appendix A for more details. Let =>■ denote 
convergence in the Skorokhod space 2?([0, 1], K) of real valued right-continuous functions with left 
limits, endowed with the J\ topology, cf. [33]. 

Theorem 4.1. Let Assumption 1 hold and assume that the function a is continuous and (5) holds 
for some q > 2a. 

(i) If 1 < a < 2 and E[Zq] = 0, then a^Sn converges weakly in the space 2}([0, 1),R) endowed 
with Skorokhod's J\ topology to an a-stable Levy process with skewness 2/3—1. 

Let t p = t(o~ p ) be the Hermite rank of the function o~ p . 
(ii) If p < a < 2p and 1 — r p (l — H) < p/a, then 

a-P(S Pin -nE{\Y Q \P})^L a/p , (18) 
where L a / p is a totally skewed to the right a/p-stable Levy process, 
(in) If p < a < 2p and 1 — t p (1 — H) > p/a, then 

n p n pl {Sp, n -nE[|YbrJ) => ; Rt p .h ■ (19) 

T p . 

(iv) If p > a, then a~ p S p . n => L a / p) where L a / p is a positive a/p-stable Levy process. 

Note that there is no effect of leverage. The situation will be different for the sample covariances. 
The fact that when the marginal distribution has infinite mean, long memory does not play any 
role and only a stable limit can arise was observed in a different context by [12]. 

5. Sample covariances 

In order to explain more clearly the nature of the results and the problems that arise, we start 
by considering the sample covariances of the sequence {Yi}, without assuming that E[Z<j] = 0. 
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For notational simplicity, assume that we observe a sample of length n + h. Assume that a > 1. 
Let Y„ = n _1 Y^j=i Yj denote the sample mean, m = E[Z ], fiy = E[Yb] = mE[<7 ] and define the 
sample covariances by 



1 ™ 

(s) = - V)(Yi - r„)(r 4+s -%), o<s<h, 



n 

i=l 



For simplicity, we have defined all the sample covariances as sums with the same range of indices 
1, . . . , n. This obviously does not affect the asymptotic theory. For s = 0, . . . , h, define furthermore 



1 " 

Cn{s) = -Y j Y i Y i+s . 

Then, defining i(s) = cov(y j Y s ), we have, for s = 0, . . . , h, 

7«0) - 7(«) = C n (s) - E[Y Q Y S ] + & - F n 2 + P (l/n) . 

Under the assumptions of Theorem 4.1, Y% — = Op(a n ). This term never contributes to the 
limit. Consider now C n (s). Recall that Ti is the sigma-field generated by (r]j,Zj), j < i and 
define 

Y _ E[X i+s I Fj-l] _ -1 

var(E[X i+s | F^}) ^ 1 + 

with <; 2 = X)j= s +i c f ■ Let if be the function defined on M 2 by 



K(x,x) = E[Z S ]E 



Z a(x)a 

c j 1 ls-j + Ss x 



E[Y Y S ] . (20) 



Then, for each i > 0, it holds that 

E[F i y i+5 1 Ji_i] - e[y y 5 ] = K(x l ,x hS ) . 

We see that if m = K[Z S ] = 0, then the function K is identically vanishing. We next write 

1 " 1 " 11 

C n (s) - E[Y Y S ] = - YiY^+s - E[YiY i+8 \ Fi-x]} + - V K (X l7 X i>s ) = -M n . s + -T n „ 

71 ± ^ 71 1 * 71 71 



n i — ' n 

i=l i=l 



The point process convergence results of the previous section will allow to prove that 6~ 1 M ns has 
a stable limit. If m = E[Z] = 0, then this will be the limit of 6~ 1 (C„(s) — E[y Xs]), regardless 
of the presence of leverage. We can thus state a first result. Let -4- denote weak convergence of 
sequences of finite dimensional random vectors. 

Theorem 5.1. Assume that a G (1, 2) and E[Zq] = 0. Under the assumptions of Propositions 3.2, 
3.3 or 3.4, 

nb-\%{l) - 7 (1), ■ • • ,%(h) - 1(h)) 4 (A, . . . ,C h ) , 
where C\, . . . ,£h are independent a-stable random variables. 
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This result was obtained by [14] in the (LM)SV case for the function <j{x) — exp(x) and under 
implicit conditions that rule out long memory. 

We continue the discussion under the assumption that m/0. Then the term T„ iS is the partial 
sum of a sequence which is a function of a bivariate Gaussian sequence. It can be treated by 
applying the results of [2]. Its rate of convergence and limiting distribution will depend on the 
Hermite rank of the function K with respect to the bivariate Gaussian vector (Xq, Xq^ s ), which is 
fully characterized by the covariance between Xq and -Xo,s> 

oo 

cov(X a , X 0jS ) = ? s -1 ^2 CjCj+s = s^Ps . 
i=i 

LMSV case Since in this context the noise sequence {Zi} and the volatility sequence {a{\ are 
independent, we compute easily that 

K(x, y) = m 2 a(x)E[a(x s C + c s i ]a + <; s y)] - m 2 E[a(X a )a(X s )} , 

where x 2 = X^=i c j an d C 1S a standard Gaussian random variable, independent of ?/o- Thus, the 
Hermite rank of the function K depends only on the function a (but is not necessarily equal to 
the Hermite rank of a). 

Case of leverage In that case, the dependence between 770 and Zq comes into play. We now have 

K(x, y) = ma(x)E[a(KsC + c s ??o + <; s y)Z ] - mE[a(X )a(X s )Z ] , 

and now the Hermite rank of K depends also on Zq. Different situations can occur. We give two 
examples. 

Example 3. Consider the case a(x) — exp(x). Then 

E[Y Y S I JU] = E[Z Z s cMXo)exp(X s ) | 7U] 



mE[Z exp(c s ?7o)]E 



J =1 



exp ^2 c jVs-j exp [X + ? s Xo, s J 



Denote m = E[Zq exp(c s ?7o)] and note that E exp ^X^=i CjVs-j^j = exp Thus 
K(x,y) = mm exp (x£/2) jexp(a; + Q,y) -E exp ( X + <; s X , s ^j | . 

If E[Zo] = or E[Zq exp (c s ?7o)] = 0, then the function K is identically vanishing and T„ ;S = 0. 
Otherwise, the Hermite rank of K with respect to (X Q ,X 0s ) is 1. Thus, applying [2, Theorem 6] 
(in the one-dimensional case) yields that n~ 1 p n 1 ^ 2 T n s converges weakly to a zero mean Gaussian 
distribution. The rate of convergence is the same as in the LMSV case but the asymptotic variance 
is different unless E[Z n cxp(c s rj )} = E[Z ]E[exp(c s ?7o)]. 
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Example 4. Consider u{x) = x 2 . Denote Xi, s — x^ 1 J2j=i c jVi+s-j- Then 

E[Y Y S | F-i] = E[Z Z s X 2 (x s X , s + ^ S X Q . S + c sVo ) 2 \ F-i] 

= mX 2 {x 2 s m + c s E[Z Q n 2 } + S s m(X , s ) 2 + 2ft,c.E[Z Jjd] Xo,*} ■ 

Thus 

K(x, y) = s s m 2 {x 2 y 2 - E[X 2 (X 2 J) + 2^c s mE[Z oVo }{x 2 y - E[X 2 X , S }} 
+ (x 2 m 2 + c s mE[Z n 2 ]){x 2 - 1) 

and it can be verified that the Hermite rank of K with respect to (Xq,Xq sS> ) is 1, except if 
E[Zo?7o] = 0, which holds in the LMSV case. Thus we see that the rate of convergence of T n , s 
depends on the presence or absence of leverage. See Example 6 for details. 

Let us now introduce the notations that will be used to deal with sample covariances of powers. 
For p > define m p = E[|Zo| p ]. If p G (a, 2a) and Assumption (1) holds, m p is finite and 
E[|Zo| 2p ] = oo. Moreover, under the assumptions of Lemma 2.1 or 2.2, for s > 0, E[|YbV,| p ] < oo 
and E[|YbV s | 2p ] = oo for p £ (a/2, a). Thus the autocovariance 7 P (s) = cov(|Yb| p , |V S | P ) is well 
defined. Furthermore, define Y p<n = n^ 1 X)"=i an d 

n 

>."( S ) = - £(W " Y p ,n)(\Y l+s \ P - Yp,n) ■ 
i=l 

Define the functions K* (LMSV case) and K^ s (case with leverage) by 

K; iS {x,y) = m y(x)E[aP(x s C + c sVQ + c s y)] - m 2 p E[a p {X Q )a p {X s )] , (21) 
Kl s {x,y) =m p a p (x)E[a p (x s C + c sVo + ^y)\Z n - m p E[a"(X )a p (X s )\Z \ p ] . (22) 

5.1. Convergence of the sample covariance of powers: LMSV case 

Theorem 5.2. Let Assumption 1 hold and assume that the sequences {rji} and {Zi} are indepen- 
dent. Let the function a be continuous and satisfy (5) with q > 4a. For a fixed integer s > 1, let 
r*(s) be the Hermite rank of the bivariate function K* s defined by (21 ), with respect to a bivariate 
Gaussian vector with standard marginal distributions and correlation c~ 7 S . 

• If p < a < 2p and 1 — r*(s)(l — H) < p/a, then 

nb~ p %, n (s) - 7 P (s) A C s , 

where C s is a a/p-stable random variables. 

• If p < a < 2p and 1 - r*(s)(l — H) > p/a, then 



Pn Tp[S)/2 {%,n{s) - 7p(s)) 4- G* , 



where the random variable G* is Gaussian if t*(s) = 1. 
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For different values s = 1, . . . , h, the Hcrmite ranks r*(s) of the functions K* s may be different. 
Therefore, in order to consider the joint autovovariances at lags s = 1, . . . , h, we define 



t* =min{r*(l),...,r*(ft)} . 



Corollary 5.2. Under the assumptions of Theorem 5.2, 

• // 1 — t*(1 — H) < p/a, then 

n K p {%A l ) - 7 P (i), ■ • -,%A h ) - 7 P (M) (A, ■ • •>£/») i 

where Ci, ■ ■ ■ ,C P are independent a/p-stable random variables. 

• If 1 — t*(1 - H ) > p/a, then 

Pn T;/ \%Al) - 7p(1). • • • . %,n(h) ~ 7p W) 4 (GJ , . . . , , 

where G* = G* if t*{s) = t* anrf G* = otherwise. 

We see that the joint limiting vector (G 1; . . . , G^) may have certain zero components if there exist 
indices s such that r*(s) > r*. However, for standard choices of the function a, the Hermite 
rank t*(s) does not depend on s. For instance, for cr(x) = cxp(a;), t*(s) = 1 for all s, and for 
cr(x) = x 2 , t*(s) = 2 for all s. 

5.2. Convergence of sample covariance of powers: case of leverage 

Theorem 5.3. Let the assumptions of Proposition 3.3 or 3.4 hold and assume that (5) holds 
for some q > 4a. Let rj(s) be the Hermite rank of the bivariate function s defined by (22), 
with respect to a bivariate Gaussian vector with standard marginal distributions and correlation 

• If p < a < 2p and 1 — rJ(s)(l — H) < p/a, then 

n K P {lpA l ) -JpO-),---,%,n(h) -l P {h)) ^ £ s , 

where C s is a a/p-stable random variable. 

• If p < a < 2p and 1 — rf(s)(l — H) > p/a, then 

p n Tp{s)/2 {%A l ) -7p(1), •■-,>,«(» -7p0)) G\ , 
where the random vector Gj is Gaussian if rj(s) = 1. 
Again, as in the previous case, in order to formulate the multivariate result, we define further 

T l =min{T 3 J(l),...,Tj(>)} . 
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Corollary 5.3. Under the assumptions of Theorem 5.3, 

• If 1 — t£(1 - H) < p/a, then 

n K P (%,ni l ) - 7 P (1), ■ • -,%,n(h) - l P {h)) A (£i, ...,C h ) , 
where C\, . . . , C p are independent a/p-stable random variables. 

• Ifl -rt(l - H) > p/a, then 

pr pV2 (>, n (l) - 7 P (1), • ■ -,%,n(h) - l P {h)) 4 (G{, . . . , G{) , 

where G\ = G\ if (s) = and G\ — otherwise. 

The main difference between Theorems 5.2 and 5.3 (or, Corollaries 5.2 and 5.3) is the Hcrmite 
rank considered. Under the conditions that ensure convergence to a stable limit, the rates of 
convergence and the limits are the same in both theorems. Otherwise, the rates and the limits 
may be different. 

Example 5. Consider the case a{x) = exp(a;). For all s > 1 we have = r^(s) = 1. Thus, 
under the assumptions of Theorem 5.3, we have: 

• If H < p/a, then nb~ {7 P ,n(s) — 7p(s)} converges weakly to a stable law. 

• If H > p/a, then p n {jp,n(s) — lp(s)} converges weakly to a zero mean Gaussian 
distribution. 

The dichotomy is the same as in the LMSV case, but the variance of the limiting distribution in 
the case H > p/a is different except if K[Zq cxp(c s ?7o)] = E[Zo]E[cxp(c s ryo)]. 

Example 6. Consider the case a(x) = x 2 and p = 1. Assume that E[t^i|Zi|] ^ 0. Then for each 
s > 1, t\ = rl(s) = 1 whereas r* = r*(s) = 2, thus the dichotomy is not the same as in the LMSV 
case and the rate of convergence differs in the case H > 1/a. 

• If H < l/a, then ^^ 1 {7n.i(s) — 7i(s)} converges weakly to a stable law. 

— 1/2 

• If if > l/a, then p n {j n ,i(s) ~ 7i(s)} converges weakly to a zero mean Gaussian 
distribution. 

If we assume now that E[7yi|Zi|] =0, then t\ = t* =2. Thus the dichotomy is the same as in the 
LMSV case, but the limiting distribution in the non stable case can be different from the one in 
the LMSV case. 

• If 2H — 1 < l/a, then n6~ 1 {7i.„(s) — 7i(s)} converges weakly to a stable law. 

• If 2H — 1 > l/a, then Pn xYl,n{ 8 ) — 7i( s )} converges weakly to a zero mean non Gaussian 
distribution. 

If moreover E[i?2 ) l-^i I] = 0, then for each s, the functions K* s and s are equal, and thus 
the limiting distribution is the same as in the LMSV case. 
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6. Proofs 

Lemma 6.1. Let Z be a nonnegative random variable with a regularly varying right tail with index 
—a, a > 0. Let g be a bounded function on [0, oo) such that \im x ^ +oc g(x) = c g 6 (0, oo). Then 
Zg(Z) is tail equivalent to Z : 

nZg(Z)> X ) = 
P(Z > x) 9 



x— > + oo 



Proof. Fix some e > and let xq be large enough so that \g(x) — c g \/c g < e for all x > x . The 
function g is bounded, thus zg{z) > x implies that z > x/||g||oo and if x > ^olMioo, we have 

V(Zg{Z) >x)= ¥(Zg(Z) >x,Z> z/||s||oo) 

< P(Zc 9 (l + e)>x,Z> x/WgWoo) < F(Zc g (l + e) > x) . 

This yields the upper bound: 

nZgjZ) > x) nZc(l + e)>x) 

nz> X ) -^fi p(z>x) = c * (1 + e) • 

Conversely, we have 

P(Z<?(Z) > x) = P(Zg(Z) >x,Z> x/WgWov) > F(Zc g (l -e)>x,Z> x/WgW*,) 
Z > xmax i — -, — — i — V I = P I Z > 



where the last equality comes from the fact that (1 — e)c g < c g — lhn s _ i . +00 g(z) < |j.g||oo- Thus 

, c nZg{Z)>x) P(Zc g (l - e) > x) 
limmt — —7 r — > hmsup = c~ 1 — e . 

x^+oo ¥{Z>x) ~ x ^+£ ¥{Z>x) gK ' 

Since e is arbitrary, we obtain the desired limit. 

Lemma 6.2. Let Z be a nonnegative random variable with a regularly varying right tail with index 
—a, a > 0. For each e > 0, there exists a constant C , such that for all x > 1 and all y > 0, 

P{yZ > g> < C(y V . (23) 

P(Z > x) ~ yy ' v ' 

Proof. If y < 1, then P(yZ > x) < P(Z > x) so the requested bound holds trivially with C = 1. 
Assume now that y > 1. Then, by Markov's inequality, 

P(yZ > z) = P(Z > x) + P(Zl {z < x} > x/y) < X - a - e y a+e E[Z a+ *l {z < x} ] . (24) 

Next, by [17, Theorem VIII.9.2] or [5, Theorem 8.1.2], 

E[Z a +n {z<x} ] a 
lim ; =-r= t- = — ■ 

x^+oo x«+ e P(Z > x) e 
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Moreover, the function x P(Z > x) is decreasing on [0, oo), hence bounded away from zero on 
compact sets of [0, oo). Thus, there exists a constant C such that for all x > 1, 

mz a+e ii Z<x x\ 

nz> X ) - • (25) 



Plugging (25) into (24) yields, for all x, y> 1, 

l + Cy 



(yZ>x) _ 1 t ^ <a+e 



V(Z > x) 
This concludes the proof of (23). 

Proof of Lemma 2.1. Under the assumption of independence between the sequences {Zi} and 
{r)i}, as already mentioned, Yq is tail equivalent to Z and YqY^ is tail equivalent to Z Zi for all 
h. The properties (11), (12), (13) are straightforward. We need to prove (14) and (15). Since Zo 
is independent of o~ 3 and Z^ by conditioning, we have 



nP(|Yo| > a n x, \Y Y 3 \ > b n x) = E 



nF\ z \ \ — V 



dnX b n X 



CO <jQ(Tj\Zj\ 

with F\z\ the distribution function of \Z$\. Since a n /b n — > 0, for any y > 0, it holds that 
lim„^ +00 nF\ z \{b n y) = 0. Thus, 

- / a„x b„x \ - ( br,x \ 
nF _JL_ v ^— < nF\ Z \ Vl H ' a ' S - 

Moreover, by Lemma 6.2 and the definition of a n , for any e > there exists a constant C such 
that 

nF\z\ (— V < nF lz{ < Cx~^a^ . 

By assumption, (5) holds for some q > a. Thus, choosing e small enough allows to apply the 
bounded convergence theorem and this proves (14). Next, to prove (15), note that \Yi\ A \Yj\ < 
(o~i V <7j ) ( | Zi | A \Zj\). Thus, applying Lemma 6.2, we have 

P(|F Y| > x, \Y Y 3 \ >x)= ¥(\Zo\ao(<Ti\Zi\ A a 3 \Z,\) > x) 

< CP(\Z \ > x)E[<j^(<j 1 V^-r+'MdZil A \Z 3 \) a +^ . 

The expectation E[crQ +e (<7i V fTj) Q+c ] is finite for e small enough, since Assumption (5) holds with 
q > 2a. Since P(|Z | > x) = o(F(\ZiZ 2 \ > x)), this yields (15) in the LMSV case. 

Proof of Lemma 2.2. It suffices to prove the lemma when the random variables Zi are non- 
negative. Under the assumption of the Lemma, exp(c^77o)^o is tail equivalent to Zq. Thus, by 
the Corollary in [16, p. 245], Z exp(ci l rio)Zf l is regularly varying with index a and tail equivalent 
to Z Q Z h . Since E[Z$] = oo, it also holds that P(Z > x) = o(P(exp(c /i 77 )Z Z 1 > x)), cf. [11, 
Equation (3.5)]. 
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Define Xh = YlkLi kjth c kVh-k- Then Xh is independent of Zq, j]q and Zh- Since Y^Yh = 
exp(Xo+X/j)Zo exp(c/j?7o)^?i, we can apply Breiman's Lemma to obtain that YqY^ is tail equivalent 
to Zoexp(cf l r]o)Zf l , hence to ZqZi. Thus (13) and (11) hold with 

d+ (h) = (i E[eMa{Xo+ * h))] , d_ W = (i-3) E[cxp(a(Xo + ^ ))] , 

M E[exp(o(X +Xi))] ^E[exp(a(X +*i))] 

where /? is the skewness parameter of .Zo exp(c/ t 7yo)^ft,- 

We now prove (15). For fixed i, j such that < i < j, define 

(oo \ / oo 

X] ^i-fc > ^t,3 = <r(Xi,j) = exp ^ Cfe^-fc 
fc=l,fc#i / \fc=l,fc^j,i-i 

Denote Zq*^ = Zoexp(cfc?7o) an d Vi = cxp(cj_i?7i). Then 

P(F ^ > x,Y Y 3 > x) = F(a Q a i Z^ Z l > x , a Q aijZ^' exp(cj-iT}i)Zj > x) 
< P( < 7 (a ?; V d-yX^ + Z^)(Z A V^Z,) > x) . 

Now, (Zi A ViZj) is independent of <7 (<3i V &i,j)(Z^ + Z^), which is tail equivalent to Zq by 
assumption and Breiman's Lemma. Thus, in order to prove (15), we only need to show that for 
some S > a, E[(Z A ViZj) 5 ] < oo. This is true. Indeed, since E[V^ 9 ] < oo for all q > 1, we can 
apply Holder's inequality with q arbitrarily close to 1. This yields for p^ 1 + q^ 1 = 1, 

E[(Zi A ViZj) 6 } < E[(l V Vi) s (Zi A Zj) 6 ] < E 1 / p [(l V V^ 8 ] E 1 /"[(Z l A Z 3 ) q& ] . 

The tail index of (Zi AZj) is 2a, and thus E 1 /«[(Z. t A^)'*] < oo for any g and 5 such that q5 < 2a. 
Thus E[(Zi A ViZ j )' 5 ] < oo for any 5 £ (a, 2a) and (15) holds. The proof of (14) is similar. 

Proof of Lemma 2.3. We omit the proof of the regular variation and the tail equivalence be- 
tween YoYh and ZqZi which is a straightforward consequence of the assumption. We prove (15). 
Using the notation of the proof of Lemma 2.2, by the subadditivity property of a , we have, for 
j > i > 0, and for some constant C , 

p(y y i > x, y y j > x) 

= V(a a(Xi + Cii] )Z Zi > x , a Q a(Xij + Cji] a + <■ , ,//,</,,/, > x} 

< V(Ca \Z \{o-{Xi) + a{cir)o)}{a(X itj ) + a(c jVQ ) + «7(c,-_<»fc-)}(|Z<| A \Zj\) > x) 

< V(Ca \Z Q \a(X t )a{X^)(\Z t \ A |Z,-|) > x) +¥(Ca \Z \a(X i )a(c j r ]0 )(\Z i \ A \Z 3 \) > x) 

+ F(Ca \Z \a(X i )a{c j . i r ]i ){\Z i \ A \Zj\) > x) +F{Ca \Z \a{c i r ]0 )a(X itj )(\Z i \ A \Zj\) > x) 
+ F(Ca \Z \o-{c i r) )a{c j r ] o)(\Z i \ A \Z S \) > x) +F(Ca \Z \a(c i r]o)a(c j - i r] i )(\Z i \ A \Z,\) > x) . 

Now, under the assumptions of the Lemma, each of the last six probabilities can be expressed as 
P(ZU > x), where Z is tail equivalent to Zq and U is independent of Z and E[|[/|' 3 ] < oo for some 



20 



Rafal Kulik and Philippe Soulier 



q > a. Thus, by Breiman's Lemma, ZU is also tail equivalent to Zo, and thus ¥(YoYi > x, Y^Yj > 
x) = O{P{\Z \ > x)) = o(F(\Y Yi\ > x)), which proves (15). 

6.1. Proof of Propositions 3.1, 3.2, 3.3 and 3.4 

We omit some details of the proof, since it is a slight modification of the proof of Theorems 3.1 
and 3.2 in [14], adapted to a general stochastic volatility with possible leverage and long memory. 
Note that the proof of [14, Theorem 3.2] refers to the proof of Theorem 2.4 in [9]. The latter proof 
uses condition (2.6) in [9], which rules out long memory. 

The proof is in two steps. In the first step we consider an m-dependent approximation X^ m ^ 
of the Gaussian process and prove point-process convergence for the corresponding stochastic 
volatility process Y^ m > for each fixed m. The second step naturally consists in proving that the 
limits for the m-dependent approximations converge when m tends to infinity, and that this limit 
is indeed the limit of the original sequence. 



First step Let x\ m ' = Y^k=\ c kVi-k, Y^' = a{x\ m ')Zi and define accordingly Y„™ . Note that 



the tail properties of the process } are the same as those of the process {Yi}, since the latter 

are proved without any particular assumptions on the coefficients Cj of the expansion (2) apart 
from square summability. In order to prove the desired point process convergence, as in the proof 
of [14, Theorem 3.1], we must check the following two conditions (which are Equations (3.3) and 



where v m is the mean measure of the limiting point process and (27) must hold for any continuous 
bounded function g, compactly supported on [0, 1] x [— oo, oo] h \ {0}. 

The convergence (26) is a straightforward consequence of the joint regular variation and the 
asymptotic independence properties (14), (15) of Yo, YqYi, . . . , Y^Yh- Let us now prove (27). Note 
first that, because of asymptotic independence, for any fixed i, 



(3.4) in [14]): 



F(Yi m 1 ) G-)A^ 



m 5 



(26) 



[n/k] 




(27) 



fim nE[ 5 (Y^) 5 (Y^)]=0. 



Next, by m-dependence, for each k, as n — > +oo, we have 



[n/k] 



[n/k] 



II 



£ E[ S ( Y W) S (Y«)] 



n 



£ E[g{Y^)]E[g(Y^)] 



i—2-\-m-\-h 



i—2-\-m-\-h 
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This yields (27). Thus, we obtain that 



h oo 



= 1 k=l 



where ^(mjj, . . . , SfcLi arc independent Poisson processes with respective mean 

measures 

A , m (dz) - a {^x""- 1 1(0,00) (a:) + (1 - /3m)(-a;)" a " 1 l(-oc,o)(^)} dx , (28) 
A s , m (dx) = a{^ n) ( S )^ Q - 1 l (0 , oo) (x)+dL m) ( S )(^)- Q - 1 l ( _ OOi0) (x)}d a; , (29) 

where d+ (s) and d_ ^(s) depend on the process considered and /3 m = /3E[er a (V( m ))]/E[<7 a (V)]. 
Second step We must now prove that 

iV m iV (30) 

as m — > +oo and that for all 77 > 0, 

lim limsupP^iVn, AK m) ) > 77) = . (31) 

TO->+O0 „_j. + oo 

where g is the metric inducing the vague topology. Cf. (3.13) and (3.14) in [14]. To prove (30), it 
suffices to prove that 

lim p m =j3, (32) 

m— s-+oo 



lim d^fs) = d+(s) , lim d { " l) (s) = d_(a) . (33) 
To prove (31), as in the proof of [14, Theorem 3.3], it suffices to show that for all e > 0, 

lim lim sup nFfa-^Yo - V M | > e) = , (34) 

m-> + oo „_^ +00 

lim lim sup nP (b~ 1 1 Y a Y s - Y~ (m) Y^ \ > e) = . (35) 

ro->+oo „^ +00 V / 

If (5) holds for some q > a and if a is continuous, then (32) holds by bounded convergence, in both 
the LMSV case and the case of leverage. We now prove (34). Since Y and Z are tail equivalent, 
by Breiman's Lemma, we have 

lim S upnP(a- 1 |y - V (m) | > <0 < Ce~ a E[\a(X^) - a(X )\ a ] . 

n— >+oo 

Continuity of a, Assumption (5) with q > a and the bounded convergence theorem imply that 
lim m ^ +00 E[|cr(^ m) ) - a{X Q )\ a ] = 0. This proves (34) in both the LMSV case and the case of 
leverage. We now split the proof of (33) and (35) between the LMSV and leverage cases. 
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LMSV case. In this case, we have 



+ w-a+w E[aa{Xo)aa{Xs)] , <t_W-*-W ______ 



For s = 1, . . . , h, define 

W m , a = a(X^ m) )a(X^) - a(X 1 )a(X 1+s ) . 

Continuity of a implies that W m ^ s — > as m — > +oo. Under the Gaussian assumption, X' m ) = 
w m X for some it TO € (0, 1), thus if (5) holds for some q' > a, then it also holds that 

supE[(T 9 '(X( m ))] < 

m>l 

hence W m converges to in L q for any q < q' . Likewise, since assumption (5) holds for some 
q' > 2a, W m ,s converges to in L q for any q < q' . Since \ W m \ and |W miS | converge to in L a , we 
obtain that d^ n '(s) and d^ n '(s) converge to the required limits. We now prove (35). Since ZqZ s 
is tail equivalent to YqYi, by another application of Brciman's Lemma, we obtain, for s = 1, . . . , h 
and e > 0, 

limBupPCft- 1 !^ -^ (m) n (m) l > e) < limsupnP(6- 1 |Z Z s ||lU miS | > Ce) < C- a e- a E[\W m J a ] 

n— t-\-oo n— !- + oo 

which converges to as m — > +oo. This concludes the proof of (35) in the LMSV case. 

To prove (35) in the case of leverage, we further split the proof between the cases <t(x) — cxp(a;) 
and a subadditive. 

Case of leverage, a(x) = exp(x) Define X s = Y^Li, ^ s c oVs-j, X^ = _"lij^ s '';'/- .. and 

W m , s = | exp(X +X S )- cxp(X { m) + X^)\ . 
As previously, we see that W m , s converges to in L q for some q > a. Thus, we obtain that 

n h oo 

E 'We?) _ — ' (n ^ +oo) ' 

i=l s=0 fe=l 

where _^Li ^r* • ■ • i _^Li <^> (m h are independent Poisson processes with respective mean 

measures A Sim (dx), s = 0, . . . ,h, defined in (28)-(29) with the constants <f+ (s) and d_(s) that 
appear therein given by 

d M (s) = d+(s) H____t^ ! d W( a )- d _( a) E ^^ m)+ ^ m) ))]. 
+ E[exp(a(X +X s ))] E[exp(a(X + X,))] 

Since |Wm, s | converges to in L q , we obtain 

oo oo 

_ 6 (tH,j k m J) _ ' ( m ^ + °°) ' * = 0, . . . , ft . 

k=l ' S k=l 
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Then, for s = 1, . . . , h, we obtain, with Zq = Zq cxp(c s ?7o), for e > 0, 

limsupnP (b-^YoYs - Y {m) Y^\ > e) = limsupnP (b- l \Z Z^ s) \\W m , s \ > e) < Ce- a E[\W m>a \ a ] 

which converges to as m — > +oo. This proves (35) and concludes the proof in the case of leverage 
with o~(x) = cxp(a;). 

Case of leverage, a subadditive We have to bound 

nP(\Z Z s \\a(X Q )a(X s ) - a(X^ m) )a(X^)\ > eb n ) . 

It suffices to bound two terms 

h(n,m) = nV(\Z Z s \\a(X ) - a(^ m) )|a(xi m >) > eb n ) , 
I 2 (n,m) = n¥(\Z Z s \a{X )\a{X s ) - a{X^)\ > eb n ) . 

Recall that xi m ^ = X^ + c s rj and X s = X s + c s t]q. By subadditivity of a, we have, for some 
constant (5, 

h{n,m) <nP(\Z Z s \\<j(X ) - a{X^ m) )\a{X^) > Ceb n ) 

+ nV(\Z Q Z s \\a(X Q ) - a(X^ m) )\a(c s r, ) > Seb n ) . 

The product ZqZ s is independent of |cr(Xo) — er(XQ m ' ) )|<7(xj m ' ) ) and tail equivalent to YqYi, thus 
we obtain 

limsupnP(|Z Z s ||a(X ) - a(X { m} )\a(X^) > 5eb n ) < Ce- a E[\cr(X ) - a(X^ n} )\ a a a (X^)} . 

n— >+oo 

We have already seen that a{X^) converges to o{Xq) in L a , thus the latter expression con- 
verges to as m — > +oo. By assumption, o-(c s r]o)\ZoZ s \ is cither tail equivalent to \ZqZ s \ or 
E[a q (c s r) )\Z Z s \ q ] < oo for some q > a, and since it is independent of |cr(X ) — £r(Xg lra ')|, we 
obtain that 

lim S u.pn¥(a(c s m)\Z Z s \\cr{X Q ) - a(X { " l) )\ > eb n ) < Ce- a E[\a(X ) - <i(X { m) )\ a } , 

n— >+oo 

where C = in the latter case. In both cases, this yields 

lim limsupnP(CT(c s 7 ?0 )|Z Z s ||cr(Xo) - cr{X { m) )\ > eb n ) = . 

Thus we have obtained that lim„ i _ i . +00 limsup Tl __j. +00 Ii(n,m) = 0. 

For the term li(n, m) we use assumption (17) with x = c s rj 0l y = X s and z = jtj; m \ Thus 

h{n,m) < n¥{\Z Q Z s \(a{c s r lQ ) V l)W m>s > eb n ) , 
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with 

W m , s = a(X ){(a(X s ) V 1) + {a{X^) V 1)}\X S - X™\ . 

Note that Wm tS is independent of \ZqZ s\(&(csVo) V 1) and W m s converges to when m — > +oo in 
L q for some q > a. Since \ZoZ s \a(c s rjo) is tail equivalent to | Yo^i | or has a finite moment of order 
q' for some q' > a, we have 

limsu P nP(|Z Z s |(a(c s 77o) V l)W m , s > eb n ) < CE[W% >S ] , 

n— ^+oo 

where the constant C can be zero in the latter case. In both cases, we conclude 
lim limsup nF(\Z Z s \(a(c s r] ) V l)W m . s > eb n ) = . 

m->+oo n^ + oo 

6.2. Proof of Theorem 4.1 

We start by studying S Pi „ . Write 

[nt] [nt] [nt] 

E m\ p - n\Y \ p }) = E - E o^n^-i]) + E (^wi^-i] - E [i y °H) 

l — l i—l 6=1 

=: M n {t) + Rn{t) . 

Note that E^l^-i] = E[|Z | p ]cr p pQ) is a function of X l and docs not depend on Z*. Then, 
by [2, Theorem 6], for r p (l — if) < 1/2 we have 

n^p-^Rn 4 ^ f /' 11 J i^ Tp , g . (36) 



T, 



I'- 



If r p (l — H) > 1/2 then by [2, Theorem 4], we obtain 

l ' 2 R n 4 ? E[|ZoHB , (37) 



where B is the standard Brownian motion and <, 2 = var(cr p (Xo)) + 2 5^°^ cov(cr p (Xo), er p (X;)). 
We will show that under the assumptions of Theorem 4.1 we have, 

< P M„ 4 L a/p . (38) 

The convergences (36), (37) and (38) conclude the proof of the theorem. We now prove (38). The 
proof is very similar to the proof of the convergence of the partial sum of an i.i.d. sequence in the 
domain of attraction of a stable law to a Levy stable process. The differences arc some additional 
technicalities. See e.g. [27, Proof of Theorem 7.1] for more details. For < e < 1, decompose it 
further as 

[nt] 

Mn{t) = E {\Yi\ P l ml<ean} ~ E [lYtfluY^ea^i-l] } 
i=l 
[nt] 

+ E {\Yi\ p l m>ean} - E [IK^I^i^I^] } =: Af« (t) + MW (f) . 
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The term M„ (■) is treated using the point process convergence. Since for any e > 0, the 
summation functional is almost surely continuous from the set of Radon measures on [0, 1] x [e, oo) 
onto 2?([0, 1], M) with respect to the distribution of the Poisson point process with mean measure 
vq (see e.g. [27, p. 215]), from Proposition 3.1 we conclude 



i=l t fc <(.) 
Taking expectation in (39) we obtain 



lim [nt]a-*>E[\Y \n {lYll>ean} ] =t f M P A (d 

J{x:\x\>e} 



(39) 



(40) 



uniformly with respect to t G [0, 1] since it is a sequence of increasing functions with a continuous 
limit. Furthermore, we claim that 



[nt] 



i=l 



o, 



(41) 



uniformly in t 6 [0,1]. We use the variance inequality (48) to bound the variance of the last 
expression by 

a r ; 2p [ni] 2 p [nt ]var(E[|r 1 |fl { | yi | >ea „ } |J- ]) < a^nt} 2 P[nt] E [(EHY^l^^To}) 21 
If p < a < 2p, by Karamata's Theorem (see [27, p. 25]) and Potter's bound, 

E[ 

Since by assumption E[cr 2a+2c (JTo)] < oo for some e > 0, for each t, we have 



[^l^l^^i^}] < Cn- 1 aP n Fz{ ^ n/(J [ x)) < Cn^ala^ix) 



[nt] 



var a-*>J2{^[\M P M\Y a \>ea n} ] - E [\Y\ P l {m>ean} \Fi-i] } 



2H-2+e,2H-t 



< Cn~ 2 [nt] 2 p [nt] < Cn ZH - 2+t t 



(42) 



where the last bound is obtained for some e > by Potter's bound. This proves convergence of 
finite dimensional distribution to and tightness in £>([0, 1],K). As in [27, p. 216], we now argue 
that (39), (40) and (41) imply that 



a n m n =*■ ^ a / p 



(43) 



and it also holds that L^j p =§■ L a / p as e —t 0. Therefore, to show (38) is suffices to show the 
negligibility of a~ p M n . By Doob's martingale inequality we evaluate 

[nt] 

sup ^X){l y il P1 {WI<"-}- E O y «l Pl {WI<«*»}l- 77 *-i]} 
te[o,i] 



< Cna- 2p E 



(\Yi\ p l {lYll<ean} - E ^Y x \v\ { \ YA<ean} \Tz\y 



<4Cna-^E[|y 1 | 2 Pl { | yi | <ea „ } ] 
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Recall that a < 2p. By Karamata's theorem (see [27, p. 25]), 

E [lY^n^y^y] ~ ^-(ea n )^F Y (ea n ) ~ ^-^"^n" 1 . (44) 

Applying this and letting £->0wc conclude that a~ p Mn^ is uniformly negligible in L 2 and so in 
probability, and thus we conclude that a~ p M n => L a / p . 

For p > a, E[|y"o| p ] = oo. In that case it is well known (see e.g. [13, Theorem 3.1]) that the 
convergence of a~ p S p . n to an a/p-stable Levy process follows directly from the convergence of 
the point process 5Zi=i ^Yi/a„ to a Poisson point process, and that no centering is needed. In the 
present context, this entirely dispenses with the conditioning argument and the long memory part 
does not appear. Therefore convergence to stable Levy process always holds. 

As for the sum S n , since E[Yo] = E[Zo] = 0, the long memory part R n is identically vanishing, 
thus in this case also only the stable limit arises. 

6.3. Proof of Theorem 5.2 

Let Ui — \YiY i+s \. We now write 



J2 - Era E[uf i ^_!]) + m i ^-i] - E[^D 

i—1 i—1 

n 

= M n . s + K;(X i} X itS ) = M n , s + T n . s . 

i—1 

As mentioned above, the second part is the partial sum of a sequence of a function of the bivariate 
Gaussian sequence (Xi,Xi jS ). The proof of the convergence to a stable law mimics the proof of 
Theorem 4.1. We split M„ iS between big jumps and small jumps. Write M„}s + Mn,s, with 

n 

M« = J2 (U?Mu,<b n c } - E[Ufl {u ^ bn£} I . 
t=i 

The point process convergence yields the convergence of the big jumps parts by the same argument 
as in the proof of Theorem 4.1. In order to prove the asymptotic negligibility of the small jumps 
parts, the only change that has to be made comes from the observation that M„}s is no longer 
a martingale. However, assuming for simplicity that we have (s + l)n observations Yi, we write, 
with Ui t k = U( s+1 )i- k = |Y( s+ i)i_feY( s+ i) i+s _fe|, 



Mi 



S Tl S 

{UitMUiKbne} ~ E Uf tk l {Ui < hne} | J"( 8+ l)i_fc_l } =: 



n,s 



Clearly, each k , k = 0, . . . , s, is a martingale with respect to the filtration {.Fj( s +i), 1 < i < 
therefore we can apply Doob's inequality and conclude the proof with the same arguments as 
previously. 
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6.4. Proof of Theorem 5.3 

Again, we mimic the proof of Theorem 4.1, however, some technical modifications are needed. 
We use the decomposition between small jumps and big jumps. To prove negligibility of the small 
jumps, we use the same splitting technique as in the proof of Theorem 5.2. To deal with the big 
jumps, the only adaptation needed is to obtain a bound for the quantity 



l^-i]) 2 l ■ (45) 



To show that (41) still holds in the present context, we must prove that the expectation in (45) is 
of order n~ 2 b 2p . The rest of the arguments to prove the convergence of the big jumps part remains 
unchanged. Note that E[|Yb}'s| p l{|Y' Y 3 |>e&„} l^-i] = G(Xo,Xo tS ), thus we need an estimate for 
the bivariate function 

G(x,y) = CT P (x)E[\Z Z s \ P a P (c s rio + Ss( + y)T-{\Z Z s \a(c s , 1(> +<; s (+y)>eb n }} > 

where £ is a standard Gaussian random variable, independent of Zo,r]o and Z s . We obtain this 
estimate first in the case <r(x) — exp(x) and then for subadditive functions. 
Let <j{x) = exp(x). As in the proof of point process convergence, we write 



Y a Y s = Z Z S cxp(c s ?7o) cxp (x + X 



By Lemma 2.2, ZqZ s exp(c s ?7o) is regularly varying and tail equivalent to ZqZ s . Since exp(p<^£) is 
independent of ZqZ s exp(c s ?yo) and has finite moments of all orders, we obtain that ZqZ s exp(c s r]o) exp(p? a £) 
is also tail equivalent to ZoZ s , hence to YqYi. Thus, by Karamata's Theorem and Potter's bounds, 
we obtain, for some <5 > 0, 

G(x, y) = exp(p(x + y))E[\Z Z s \ p exp(pc a r)o) cxp(p^ s C)l{\z„z s \ cx P (p Cs r, ) exp(f.O)>e&„ ex P (-»)}] 
< Cn~ x b p n exp(px) exp((p - a + S)(y V 0)) . 

Since the log-normal distribution has finite moments of all order, we obtain that E[G 2 (A , X , s )] = 
0(n~ 2 h 2 v^ which is the required bound. This concludes the proof in the case <j(x) = exp(ic). 

Let now the assumptions of Proposition 3.4 be in force. Using the subadditivity of a p , we 
obtain G(x, y) < Yn=i Ii( x > v) witn 

h{x,y) = cT p (x)E[\Z Z s \ p a p ($ s )l {l z z B \*(* s )>eK}l 
I 2 (x,y) = a p (x)E[\Z Z s \ p a p (y)l {]ZQZsW(y)>tbn} ], 
I 3 (x,y) = cr p (x)E[\Z Z s \ p a p (d s )l UZoZsltT{y)>ebn} ], 
h(x,y) = a p (x)E[\Z Z s \ p a p (y)l {lZoZs ^ 3)>ebn} } , 

where for brevity we have denoted i9 s = c s ?7o + <; s £. We now give bound E[Ij(X ,X ^ s )], j = 
1, 2, 3, 4. Since by the assumptions, \Z^Z s \a{ , d s ) is tail equivalent to \ZqZ s \, Karamata's Theorem 
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yields 

<jP(x)E[\Z Z s \P<jP(V s )l {lZoZsW{i)s)>ebn} } < Cn-^a^x) , 

and since E[<r 2p (Xo)] < oo by assumption, we obtain by integrating that E[I 2 (Xq, Xq, s )} = 
0(n~ 2 b 2 P). For I2, using again Karamata's Theorem and Potter's bound, we obtain, for some 
8 > 0, 

aV(x)E[\Z Z s \Pa*(y)l {lZoZsHy)>ebn} ) < Cn-^^Hv) V l) p ~ a+d • 

Since |Z |cr(i? s ) is tail equivalent to \Z \ and Z s is independent of Z (t(i? s ), we easily obtain a 
bound for the tail of \Z Z s \{a{d 8 ) V 1): 

F(\Z Z s \(a(tf s ) Vl)>x)< ¥{\Z Q Z s \a{d s ) > x)+F{\Z Q Z s \ > x) < CP(Z Z S > x) , 

for x large. Thus, applying Karamata's Theorem and Potter's bound to \ZqZ s \ yields, for some 
arbitrarily small 6 > 0, 

h(x,y) < C<jP(x)E[\Z Z s \n {aiy)lZoZsl>ebn} ] < Cn-X^WWv) V 1) Q+5 

and thus we conclude that E[I|(X , X 0tS )] = 0{n- 2 b 2 P). Finally, we write, 

h(x,y) < aP(x)aP(y)E[\Z Z s \P K(tf s ) V 1) 1 { \z oZb \( 

and by the same argument as for I3 we obtain that E[I 2 (Xo, Xo tS )] = Oin^b 2 ^). 

Appendix A. Gaussian long memory sequences 

For the sake of completeness, we recall in this appendix the main definitions and results 
pertaining to Hcrmite coefficients and expansions of square integrable functions with respect 
to a possibly non standard multivariate Gaussian distribution. Expansions with respect to the 
multivariate standard Gaussian distribution are easy to obtain and describe. The theory for non 
standard Gaussian vectors is more cumbersome. The main reference is [2]. 

A.l. Hermite coefficients and rank 

Let G be a function defined on R k and X = (X^\ . . . , Jfw) be a fc-dimensional centered 
Gaussian vector with covariance matrix T. The Hermite coefficients of G with respect to X are 
defined as 



J(G,X,q) 



k 

G(X)l[H qj (X^) 
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where q = (q±, , . , , <?/.) £ N k . If T is the k x k identity matrix (denoted by Ik), i.e. the components 
of X are i.i.d. standard Gaussian, then the corresponding Hermitc coefficients are denoted by 
J*(G,q). The Hermite rank of G with respect to X, is the smallest integer r such that 

J(G, X, q) = for all q such that < \q\ H h q k \ < t . 

A. 2. Variance inequalities 

Consider now a /c-dimcnsional stationary centered Gaussian process {X^, i > 0} with covariancc 
function p n (i, j) = ^[X^ Xn^} and assume either 

oo 

VI < i,i < fc , J2\Pn(i,j)\ <oo, (46) 
or that there exists H G (1/2, 1) and a function I slowly varying at infinity such that 



Pnihj) 

n-n-oo n 2H ~ 2 £(n) 

and the bjj-s are not identically zero. Denote then p n = n 2H ~ 2 l(n). Then, we have the following 



cases. 

• If (47) holds and 2r(l — H) < 1, then for any function G with Hermite rank r with respect 
to X , 

var ^n- 1 ^ G(Xj) j < Cp^ var(G(X )) . (48) 

• If (47) holds and 2r(l — H) > 1, then for any function G with Hermite rank r with respect 
to X , 

var g ( x j) j < var(G(X )) . (49) 

• If (46) holds, then (49) still holds. 

In all these cases, the constant G depends only on the Gaussian process {X^} and not on the 
function G. The bounds (48) and (49) are Equation 3.10 and 2.40 in [2], respectively. The 
bound (49) under assumption (46) is a consequence of Equation 2.18 in [2, Theorem 2]. 

A. 3. Limit theorems 

We now recall [2, Theorem 6]. Let again {X;} be a stationary sequence of fc-dimensional 
Gaussian vectors with covariancc matrix G and such that (47) holds, and let r be the Hermite 
rank of G w.r.t. Xo. If r(l — H) < 1/2, there exists a process Rg,t.h such that 

, [«•] 

— — ^(G(X i )-E[G(X )])^ J R G>T>ff . (50) 
npn f= i 
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In particular, if k = 1, then 

1 [ ™' ] T (C\ 

— 2 £{G(A,) - E[G(A )]} 4 ^Yr t .h , (51) 



where J T (G) = ¥,[G(Xi)H T (Xi)] and Rt,h is the so-called Hermite or Rosenblatt process of order 
r, defined as a r-fold stochastic integral 

exp (it(xi + • • • + x T )) - 1 -q 

i=l 

where W is an independently scattered Gaussian random measure with Lebesgue control measure 
and 

Kf(r,H) = 



RrMt) = *Ti(r, ff) [°° ... [°° expdt^ + .-.+xO)-! jj . . . ^ (dXr) , 

J-oo j— oo xi + • • • + a; T 



m _ (T(ff-l) + l)(2r(H-l) + l) 



r! {2r(2 - 2#) sin7r(iJ - ±)} r ' 
In particular, for r = 1, then the limiting process is the fractional Brownian motion, which is a 
Gaussian process, so 



1 ™ / 

£{gpq)-e[g(x )]}4n o 



J(l) 



1/2 2L-r«— V V L-V u/jj I ' H(2H — 1) 

On the other hand, if 1 - r(l - H) < 1/2, then 

1 M 

-^^{G(A 4 )-E[G(A )]}^B, (52) 
v n i=l 

where B is the standard Brownian motion and q 2 = var(G(Xo)) + 2 cov (G(Ao), G(X,)), the 

latter series being absolutely summable. 
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